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ABSTRACT 

Measurements  of  temperature,  wave  height,  and  orthogonal  water 
particle  velocity  were  made  in  May  1973  at  the  NUC  Tower  located  one 
mile  off  Mission  Beach,  California.   A  detailed  analysis  of  the  tem- 
perature field  was  made  using  digital  temperature  and  isotherm  contour 
plots.   Billow  turbulence  microstructure  appeared  to  be  the  dominant 
mechanism  with  the  exception  of  one  run  in  which  evidence  of  double- 
diffusion  microstructure  was  found.   Spatial  correlation  lengths, 
calculated  from  the  plot  of  the  covariances,  were  of  the  order  of 
130  cm  and  less  when  the  signal  was  high-pass  filtered  for  waves  of 
100  seconds  and  longer.   No  depth  dependence  was  noticed.   Both  fre- 
quency and  wavenumber  spectra  were  calculated  and  a  correspondence 
between  the  spectra  was  noted  at  the  frequency  and  wavenumber  of  the 
surface  wave-induced  particle  displacements. 

The  Thornton,  Boston,  Whittemore  model  of  wave-induced  tempera- 
ture fluctuations  was  tested  and  found  to  model  the  temperature  spec- 
tra quite  well,  especially  in  a  narrow  band  of  frequencies  associated 
with  surface  waves.   The  turbulent  temperature  spectrum,  calculated 
as  the  difference  between  the  actual  and  wave-induced  spectra,  had  a 
slope  near  -5/3  above  0.1  Hz  and  more  negative  at  lower  frequencies. 
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I.    INTRODUCTION 

A.   HISTORY 

For  the  past  three  years  research  has  been  carried 
out  at  the  United  States  Naval  Postgraduate  School  (NPS) , 
Department  of  Oceanography,  to  study  the  shallow  water 
interactions  that  occur  between  various  oceanographic 
parameters  such  as  temperature,  salinity,  surface  waves, 
internal  waves,  and  sound  speed.   These  experiments  have 
been  conducted  from  the  Naval  Undersea  Research  and  Devel- 
opment Center's  (NUC)  Oceanographic  Research  Tower  located 
approximately  one  mile  off  Mission  Beach,  California. 
Results  of  the  experiments  are  presented  in  NPS  Theses, 
NPS  technical  reports  and  publications.   This  thesis  is 
based  on  a  study  started  by  Whittemore  (1973)  and  is  an 
extension  of  the  results  obtained  by  Hinard  (1973). 

The  Whittemore  experiment  was  conducted  at  NUC  on  8 
and  9  June  1972.   Measurements  were  made  of  salinity, 
sound  velocity,  particle  velocity,  and  temperature  (at 
various  positions  along  a  horizontal  array)  for  ten  pre- 
selected depths.   He  concluded  that  the  relationship 
between  temperature  fluctuations  and  surface  wave  height 
(for  swell  type  waves)  was  of  the  form, 


T(t)  =  K(f,z) 


dT 
dz 


n(t)  , 


and  that  a  similar  relationship  was  likely  to  exist  be- 
tween temperature  fluctuations  and  internal  waves.   He 
also  found  that  turbulence  degrades  coherence  between 
temperature  fluctuations  in  sensors  displaced  spatially; 
however,  a  lack  of  temperature  sensors  prevented  him 
from  making  definitive  statements. 

Reworking  the  data  collected  during  the  June  1972  ex- 
periment, Thornton  (1974)  demonstrated  a  means  for  sepa- 
rating turbulent  and  surface  wave-induced  vertical  veloc- 
ity spectral  components  allowing  for  the  quantities  to 
be  statistically  correlated.   The  method  applied  to  the 
wave  and  velocity  measurements  allowed  the  wave  energy 
density  spectral  components  to  be  converted  to  velocity 
spectral  components  using  linear  wave  theory.   The  com- 
puted values  compared  very  well  with  the  measured  veloc- 
ity spectra  for  the  moderate  wave  conditions  encountered 
and  appeared  to  verify  the  theoretical  spectral  transfer 
function  and  assumed  linear  system. 

Using  the  same  data  and  techniques,  Thornton,  Boston, 
and  Whittemore  (1974)  devised  a  means  for  separating  tur- 
bulent and  surface  wave-induced  temperature  fluctuations. 
They  assumed  that  linear  wave  theory  could  be  used  to  de- 
scribe the  surface  waves;  that  the  temperature  isotherms 
were  horizontal;  and  that  the  mean  temperature  gradient 
was  a  constant.   In  this  way  they  developed  a  stochastic 
model  describing  the  surface  wave-induced  temperature 
fluctuations.   This  allowed  a  means  for  calculating  the 
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turbulent  temperature  spectrum  by  subtracting  out  the 
wave-induced  contributions.   The  measured  slope  of  the 
turbulent  temperature  spectra  closely  approximated  the 
Kolmogorov  -5/3  slope  over  a  decade.   This  supported  the 
result  that  the  turbulent  temperature  fluctuations  are 
statistically  independent  of  the  surface  wave-induced 
velocities . 

In  the  present  work  it  is  intended  to  refine  the 
analysis  of  the  data  and  extend  the  results  obtained  by 
Minard  (1973).   Minard  set  out  to  rectify  the  lack  of 
sensors  noted  in  Whittemore's  experiment  and  conducted 
another  experiment  at  NUC  on  15  and  16  May  1973.   Measure- 
ments of  small  scale  temperature  fluctuations  were  made 
in  19  meters  of  water  using  a  movable  array  of  seven 
thermistors,  with  a  short  (0.015  sec)  response  time;  and 
accurate  to  0.01°C  in  seawater.   A  schematic  of  the  tower 
and  instrument  mounting  is  given  in  Figure  1.   In  addi- 
tion to  temperature,  wave  heights  and  orthogonal  water 
particle  velocities  were  measured  simultaneously.   Mea- 
surements were  made  at  various  depths  with  either  a  ver- 
tical or  horizontal  orientation  of  the  temperature  array. 

Minard  found  that  in  the  presence  of  internal  waves 
temporal  scales  were  of  the  order  of  4  seconds  and  spatial 
scales  were   of  the  order  of  33  cm,  whereas  outside  the 
influence  of  internal  waves  they  were  of  the  order  of  10 
seconds  and  118  cm.   It  was  concluded  that  in  the  presence 
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Figure  1.   Schematic  of  Tower  and  Instrument  Mountinc 
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of  internal  waves,  correlation  times  were  similar  for 
both  vertical  and  horizontal  scales  and  spatial  corre- 
lation lengths  were  similar  for  vertical  and  horizontal 
scales . 

B.   OBJECTIVES 

The  main  objectives  of  this  thesis  with  relation  to 
the  measurements  made  by  Minard  in  May  1973  are  to: 

1.  Define  the  spatial  and  temporal  temperature 
field.  . 

2.  Define  the  spatial  and  temporal  temperature 
spectra  and  determine  if  an  obvious  relationship  exists 
between  them. 

3.  Define  the  temperature  gradient  field  and  its 
variation  with  time  to  enable  a  comparison  to  be  made 
between  the  formation  of  microstructure  with  and  with- 
out the  passage  of  an  internal  wave. 

4.  Test  the  stochastic  model  proposed  by  Thornton, 
Boston,  and  Whittemore  to  calculate  the  theoretical 
wave-induced  temperature  spectrum  using  these  data. 
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II.   THEORETICAL  CONSIDERATIONS 

The  creation  of  small  scale  temperature  fluctuations 
in  shallow  water  is  a  complex  process  occurring  at  all 
depths  and  scales,  depending  on  the  operative  driving 
function  prevailing  in  time  and  space.   In  the  near  sur- 
face regime  winds  over  the  surface  produce  tangential 
surface  stresses  both  from  the  interfacial  stress  and 
from  the  momentum  loss  associated  with  processes  such  as 
breaking  waves.   These  result  in  a  downwards  momentum 
transfer  with  its  attendant  downward  flux  of  heat  and 
salt  from  the  surface  to  a  depth  dependent  on  the  sta- 
bility of  the  water  column,  the  gradient  of  shears  pro- 
duced, and  the  wavelength  of  the  surface  waves.   In  this 
way  the  temperature  and  salinity  in  the  mixed  layer  be- 
come virtually  uniform;  and,  unless  the  stability  or 
Brunt-Vaisala  frequency  (defined  as 


N2  =  -  g  9P 


p  dz     ' 


(1) 


where  g  is  the  acceleration  of  gravity,  p  is  the  density 
and  ■* —  is  gradient  of  density)  is  less  than  or  equal  to 
zero,  a  thermocline  develops  due  to  the  increasing  con- 
trast between  the  properties  of  the  water  in  and  below 
the  mixed  layer.   If  there  is  substantial  surface  cooling, 
there  may  be  a  region  below  the  mixed  layer  in  which  the 
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stability  frequency  is  less  than  zero  and  convective  mo- 
tions can  develop  extending  to  considerable  depth,  parti- 
cularly in  polar  waters  (Phillips,  1969). 

The  ocean  bottom  exerts  a  retarding  force  on  any  near 
bottom  water  moving  over  it.   The  resulting  shear  zone 
again  causes  turbulent  mixing.   In  the  absence  of  dynamic 
instabilities  leading  to  a  complete  overturn,  the  bottom 
turbulence  will  produce  a  mixed  layer  similar  to  that 
formed  in  the  surface  layer  and  consequently  intensify  the 
the  thermocline  from  the  bottom. 

If  a  thermocline  is  formed,  it  will  inhibit  the  verti- 
cal extent  of  mixing  from  either  above  or  below  depending 
on  the  depth  of  the  thermocline  and  the  intensity  of  the 
gradient  within  it.   Turbulence  produces  changes  which 
tend  towards  three-dimensional  isotropy,  and  as  such  will 
be  seriously  inhibited  by  a  constraint  in  any  direction. 
In  this  case,  the  vertical  motion  is  inhibited  and  re- 
sults in  a  flow  of  energy  out  of  the  other  two  components 
of  velocity  in  an  effort  to  boost  the  vertical  energy. 
Hence,  an  increase  in  stability  as  defined  by  an  increased 
degree  of  stratification  or  density  gradient  may  well  be 
associated  with  the  degradation  of  turbulence  within  the 
water  column,  and  in  particular,  that  associated  with 
internal  waves. 

Paradoxically,  the  thermocline,  which  tends  to  inhibit 
turbulent  mixing  from  the  vertical  boundaries,  may  itself 
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be  a  source  of  increased  turbulence.   Internal  waves 
propagate  along  the  thermocline  and  can  lead  to  increased 
turbulence  through  shears  developed  between  water  masses 
either  above  or  below  the  thermocline.   In  addition, 
under  certain  conditions  the  different  rates  of  diffusion 
of  heat  and  salt  across  the  thermocline  lead  to  a  small 
scale  perturbation  in  the  temperature  field. 

For  the  purposes  of  this  thesis,  the  scales  to  which 
turbulent  temperature  fluctuations  may  be  broken  down 
fall  into  three  convenient,  though  not  necessarily  inde- 
pendent, classifications.   These  are:  micros tructure 
effects,  surface  wave  effects,  and  internal  wave  effects. 
Each  of  these  will  be  discussed  separately. 

A.    MICROSTRUCTURE  EFFECTS 

Much  evidence  is  now  available  which  shows  that  ver- 
tical mixing  in  the  interior  of  stable  fluids  occurs, 
on  occasion,  at  scales  much  smaller  than  the  vertical 
extent  of  the  water  column  (Woods  and  Wiley,  1972) . 
Although  turbulence  can  enter  at  any  boundary  of  a  fluid, 
as  long  as  vorticity  is  generated,  two  processes  have 
been  proposed  as  capable  of  producing  intense  vertical 
micros tructure  activity.   The  first  arises  from  the  differ- 
ing rates  of  diffusion  of  heat  and  salt  (double-diffusion 
or  salt-fingering).   The  second  process  is  turbulent  mix- 
ing driven  either  by  surface  winds  or  by  shear  instabili- 
ties (billow  turbulence).   It  occurs  at  the  vertical 
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boundaries  or  on  internal  waves  at  depth. 

1 .   Double- Diffusion  and  Salt-Fingering 

This  process  is  a  local  one  operating  where  weak 
mixing  has  already  increased  the  vertical  gradients  of 
temperature  and  salinity,  and  whose  vertical  scale  seldom 
exceeds  a  few  meters  (Gregg,  1973).   Considering  only  a 
two  layer  system  involving  both  heat  and  salt  gradients, 
three  distinct  situations  can  arise. 

The  first  is  an  inherently  stable  system  where 
the  upper  layer  is  both  warmer  and  less  saline  than  the 
lower  layer.   It  is  of  less  interest  because  both  compo- 
nents act  together  to  produce  a  stable  density  step 
across  the  interface. 

The  second  situation  is  an  inherently  unstable 
system  where  the  upper  layer  is  both  warmer  and  more 
saline  than  the  lower  layer.   In  a  very  few  minutes  there 
will  be  a  rapid  loss  of  heat  in  the  upper  layer  making  it 
more  dense  than  the  lower  layer,  which  is  simultaneously 
gaining  heat  and  consequently  becoming  less  dense.   The 
result  is  alternate  fingers  of  'heavier'  water  sinking 
and  'lighter'  water  rising,  hence  its  name  -  Salt-Finger- 
ing (Figure  2) .   The  pattern  of  convection  will  persist 
only  to  a  limited  depth  where  the  descending  more  saline 
water  undergoes  weak  overturning  to  form  a  well  mixed 
layer  of  temperature  and  salinity  intermediate  to  the 
parent  waters.   In  such  a  system  a  series  of  steps  will 
form  beneath  the  interface  with  steps  of  the  order  of 
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Figure  2.   SALT-FINGERING-  The  solid  lines  show  the  ini- 
tial and  subsequent  density  profiles  when  the 
more  rapid  diffusion  of  heat  has  produced  a 
potentially  unstable  situation  at  the  inter- 
face. (From  Gregg,  1973) 
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raised    in    the    lower 
part    of    the    record    illustrat- 
ing   the    disturbance    caused 
by    the    wire    and    the    probe 
moving    through    the    water    ahead 
of    the    sensor. 


340m 


340m 


330m 


Figure    3.      High  resolution   temperature  profile  of   the    310-345m 

interval    in   the  Arctic  off   Ice-Island   T-3    (From  Denner, 
1971) . 
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20-30  centimeters  and  interfaces  as  small  as  two  centi- 
meters . 

The  third  situation  is  one  in  which  the  layers 
remain  separate  but  where  convective  mixing  occurs  within 
each.   It  is  characterized  by  a  cooler,  less  saline  upper 
layer  and  a  warmer,  more  saline  lower  layer.   Again  there 
is  a  rapid  diffusion  of  heat  across  the  interface  but  now 
it  acts  to  increase  the  stability  between  the  two  layers 
with  intense  convective  mixing  within  each  of  the  layers. 
In  a  manner  similar  to  that  for  salt- finger ing  the  convec- 
tion is  self -limiting  and  extends  into  a  series  of  steps 
above  the  original  interface.   Neshyba,  Neal  and  Denner 
(1971)  found  as  many  as  40  steps  above  a  layer  of  warm 
saline  Atlantic  water  that  had  entered  the  Arctic.   Each 
is  about  3-10  meters  per  step  and  with  interfaces  of  the 
order  of  centimeters  (Figure  3) . 
2 .   Billow  Turbulence 

Billow  turbulence  is  described  as  free  shear  tur- 
bulence modified  by  a  density  gradient  and  initiated  by 
Kelvin-Helmholtz  instability  (Woods  and  Wiley,  1972) 
which  is  illustrated  in  Figure  4.   This  micros  true ture  is 
very  complex  in  form  as  opposed  to  the  rather  simple, 
symmetrical  configuration  of  the  density  steps  in  double- 
di  f f us  ion . 

It  appears  that  a  great  majority  of  the  flow  with- 
in the  thermocline  is  laminar  and  devoid  of  turbulent 
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Figure  4.    KELVIN-HELMHOLTZ  SHEAR  INSTABILITY 
The  upper  layer  has  a  higher  velocity  and  a  lower  density 
than  the  lower  layer.   The  instability,  generated  in  a 
two-layer  laboratory  tank  by  rapidly  raising  and  then 
lowering  one  end  of  the  tank,  caused  the  two  layers  to 
soon  lose  their  coherence  and  break  up  into  turbulent 
patches.   The  diagrams  are  based  on  a  study  by  S.  A. 
Thorpe  of  the  National  Institute  of  Oceanography  in 
England.    (From  Gregg,  1973) 
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motions  on  any  scale  (Woods  and  Wiley,  1972).   Occasion- 
ally, however,  shear  instability  is  triggered  within 
small  regions  by  boundary  produced  shears  or  internal 
gravity  waves.   In  this  latter  case,  as  one  of  the  wave 
packets  travels  along  a  sheet,  shear  instability  is  pro- 
duced at  the  crests  of  sufficient  height  to  cause  the 
local  Richardson's  number  of  the  billow  to  fall  below  a 
critical  value  of  0.25  (Woods  and  Wiley,  1972).   The 
Richardson  number  for  the  sheet  is  given  by 


J  = 


ghAp 
pAV2 


where  "g"  is  gravitational  acceleration,  "h"  is  the  thick- 
ness of  the  shear  zone,  "Ap/p"  is  the  fractional  density 
difference,  and  "AV  is  the  velocity  difference  across  the 
shear.   For  a  given  value  of  J  <  0.25  the  layer  will  be 
unstable  to  a  specific  range  of  wavelengths,  of  which 
L   =  7.5  x  h,  will  be  the  fastest  growing  disturbance 
(Woods ,  1968a) . 

The  condition  for  termination  of  turbulent  motions 
is  that  the  Richardson  number  again  increases  above  the 
critical  value  of  0.25  to  about  unity.   When  this  is 
reached,  the  billow  no  longer  extracts  sufficient  kinetic 
energy  from  the  shear  to  supply  the  billow  turbulence  and 
the  turbulent  motions  decay.   Woods  and  Wiley  (1972)  found 
that  in  the  sea  the  approximate  scales  were  75  cm  for  the 
billows,  500  cm  for  the  patch  of  turbulence,  and  300  sec 
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for  the  lifetime  of  the  turbulence. 

B.    SURFACE  WAVE  EFFECTS 

During  the  past  few  years  considerable  effort  has 
been  put  into  understanding  the  effects  of  wave  motion  on 
the  thermocline  but  relatively  little  has  been  done  to 
understand  the  internal  temperature  effects  such  waves 
produce.   Clearly,  temperature  fluctuations  in  the  near 
surface  region  can  be  caused  by  a  number  of  factors,  in- 
cluding surface  waves  and  their  associated  turbulence, 
internal  waves,  and  advection.   Surface  wave  effects  on 
temperature  fluctuations  are  usually  considered  only  to 
depths  of  less  than  one-half  of  one  wavelength  due  to 
the  exponential  decay  of  particle  motion  with  depth.   How- 
ever, the  energy  densities  of  surface  wave-induced  fluc- 
tuations are  in  the  same  frequency  band  and  of  the  same 
order  of  magnitude  when  they  are  present  simultaneously 
and  hence  are  very  difficult  to  separate.   Therefore  the 
spectral  separation  of  the  turbulent  and  wave-induced 
temperature  fluctuations  will  be  considered  in  greater 
detail  later  in  this  thesis. 

The  temperature  field  can  be  considered  as  a  turbulent 
field  having  an  "energy"  containing  region,  an  "inertial" 
subrange,  and  an  "energy"  dissipation  region  analogous  to 
that  hypothesized  for  velocity  fields.   In  addition,  since 
the  rates  of  thermal  and  viscous  dissipation  arc  different 
in  the  ocean,  the  smaller  scale  temperature  field  will 
tend  to  persist  longer  than  the  velocity  field. 
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1 .   Surface  Wave-induced  Temperature  Fluctuations 

Earlier,  mention  was  made  of  a  paper  (Thornton  et 
al,  1974)  setting  out  a  relationship  between  wave-induced 
temperature  fluctuations  and  the  waves.   The  major  portion 
of  what  follows  is  taken  directly  from  that  paper  and  is 
included  here,  in  a  somewhat  abbreviated  form,  for  com- 
pleteness . 

Thornton  (1974)  had  previously  shown  that  linear 
wave  theory  estimated  the  wave-induced  water  particle 
motion  very  well  under  small  amplitude  wave  conditions. 
Since  it  worked  for  particle  motion,  linear  wave  theory 
was  assumed  to  estimate  the  wave-induced  temperature  fluc- 
tuations.  This  implies  that  the  flow  is  irrotational . 
It  was  further  assumed  that  the  temperature  fluctuations 
are  small,  allowing  density  to  be  considered  constant  over 
short  distances.   Hence,  temperature  is  considered  to  be 
a  passive  scalar  with  buoyancy  effects  neglected.   Assum- 
ing temperature  is  a  conservative  property,  at  least  in 
the  body  of  the  flow,  the  conservation  of  temperature 
flux  is  given  by 


36        36  36 

-r—  +   U  .   -5—  +   V  -it— 

3 1     1  dx .  3z 

1 


(2) 


where  6(x.,z,t)  is  the  measured  temperature ,  u . (x . , z , t ) 
are  the  horizontal  velocities,  v(x. ,z,t)  is  the  vertical 
velocity,  x.  are  the  horizontal  cartesian  coordinates, 
z  is  the  vertical  cartesian  coordinate.   The  equation 
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considers  the  individual  time  change,  advection,  and  dif- 
fusion of  the  passive  scalar  temperature.   In  terms  of 
their  means,  turbulent  fluctuations  and  wave-induced 
fluctuations  the  quantities  are  defined  as 


e    =  6  +  6   +8* 
w 


u.=  u.+  (u  ).  +  u  .  ' 

11       W  1       1 


V  =  V  +  V  ' 

w 


where  the  overbar  indicates  a  time-averaged  mean,  the  sub- 
script w  indicates  contributions  from  surface  waves,  and 
the  prime  indicates  turbulent  motions. 
Assuming  that 

(1)  turbulent  and  wave-induced  motions  are  sta- 
tistically independent, 

(2)  mean  temperature  isotherms  are  horizontal, 

(3)  mean  velocities  are  zero,  and 

(4)  the  mean  vertical  temperature  gradient  is 
constant , 

allows  separation  of  the  turbulent  motions  and  Equation 
(2)  reduces  to 


80  36 

-5-^  +  (u  ).  3-^  +  v   -_  (0  +  0  )  =  0 
dt       widx.     w3z        w 


(3) 


where  the  velocities  and  the  assumption  of  irrotational 
flow  are  as  specified  in  linear  wave  theory. 

Temperature  fluctuations  caused  by  wave-induced 
velocities  are  a  priori  assumed  equal  to  zero  and 
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therefore 


dQ  89 

(u  )  .  -5—^  +  V   3-5L  =  0  . 

w  1  dx  .  w  3z 

1 


This  leaves  a  first  order,  constant  coefficient, 
homogeneous  differential  equation,  if  the  mean  tempera- 
ture gradient  is  constant,  and  therefore 


d9  r,  5 

88   V   .  (4) 


dt       9z    w 


This  is  easily  integrated  to  give 


q     /               4-\           r       38    smh    k(d  +  z),  .,  „    ..  ,_. 

6     (x.,z,t)     =     [-    -5 : — -         ]     a    cos     (k.x.-a.t)(5) 

w       1  dz    smh    kd  111 


where  V  (x. ,z,t)  is  the  wave-induced  vertical  water-par- 
w   1 

ticle  velocity  specified  by  linear  wave  theory,  'a'  is 
the  surface  wave  amplitude,  k  is  the  wave  number,  0    is 
the  radian  frequency,  d  is  the  depth,  and  z  is  the  nega- 
tive downwards  displacement  from  the  mean  water  surface. 

In  terms  of  the  spectral  transfer  function,  Hfl (0) 
(the  square  bracket  in  Equation  (5)),  and  the  sinusoidal 
surface  displacements,  r)(x.,t),  Equation  (5)  reduces  to 


)  (x.,z,t)  =  Hfl  (a)  nU.  it)  (6) 

W    1  U  X 


Equation  (6)  satisfies  Equation  (3)  exactly,  which  means 
that  moderate  swell  type  waves  do  little  or  no  mixing, 
but  only  serve  to  pump  the  thermocline  up  and  down. 
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2 .   Turbulent  Temperature  Spectrum 

In  a  manner  similar  to  that  discussed  in  Thornton, 
(1974)  the  paper  by  Thornton,  Boston,  and  Whittemore  (1974) 
discusses  the  separation  of  the  turbulent  temperature 
spectrum  from  the  wave-induced  temperature  spectrum.  For 
completeness  the  arguments  are  repeated  here. 

The  wave-induced  temperature  spectrum,  SQ  (0) ,  is 

w 
calculated  from  the  wave  spectrum,  S  (a) ,  using  the  com- 
plex transfer  function,  Hfl(a)   (given  in  Equation  (6)),  in 


the  form 


sfl  (a)  =  |h„  (a)  |  "  s    la)    . 


(7) 


w 


In  terms  of  the  total  temperature  spectrum  the 
turbulent  temperature  spectrum  is  given  by 


sQl(a)  =  sQ(a)+se  (a)-2Re{HQ (a) «sQ  (a)}  .        (8) 

w  n 


If  the  cross-spectrum  is  expressed  in  terms  of 
its  co-  and  quadrature  spectra 


S0  (a)  =  CQ  (o)  +  iQQ  (a) 

n  n  n 


then  Equation  (8)  becomes 


s01(o)  =  sQ(a)  +  se  (o)  -  2Hg(o)  Cg  (a)  (9) 


w 


which  gives  the  turbulent  temperature  spectrum  in  terms 
of  measurable  quantities. 
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Assuming  statistical  independence  between  turbu- 
lent and  wave-induced  temperature  leads  to  further  sim- 
plification, since 


SQ  (a)  =  SQ   (a) 

n       wti 


and 


HA(a)sfl   (a)  =  s   (a)  =  | h„  (a)  |   s  (a) 


'9  l"'"6 


wn 


w 


n 


This  allows  the  turbulent  temperature  spectrum  to  be 
given  simply  as  the  actual  temperature  spectrum  minus  the 
wave-induced  temperature  spectrum,  or, 


s01  (a)  =  sQ  (a)  -  |hq  (a)  |   s  (a) 


(10) 


Further  the  coherence  between  temperature  fluctua- 
tions and  the  waves  is  then  given  by 


2  se-(a)   -i 

^e     (a)  "   [1  +  irHzY] 
n 


sQ  (a) 


w 


sQ  (a)  * 


(11) 


w 


C.    INTERNAL  WAVE  EFFECTS 

Internal  waves  are  subsurface  waves  existing  between 
layers  of  different  density  or  within  layers  where  a  ver- 
tical density  gradient  exists.   They  occur  often  due  to 
a- variety  of  causes  such  as  flow  over  an  irregular  bottom, 
atmosphere  disturbances,  tidal  forces,  and  shear  flow. 
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The  theory  of  the  existence  of  internal  waves  implies 
perfect  coherence  of  plane  waves  over  space;  whereas  in 
practice  naturally  occurring  internal  waves  are  not  so 
coherent.   In  addition,  the  ocean  contains  a  number  of 
generating  sources  which  may  reinforce  or  cancel  each 
other  (LaFond,  1966) . 

In  shallow  water,  at  the  same  site  as  this  experiment, 
LaFond  (1966)  found  internal  waves  of  amplitudes  ranging 
from  2  to  23  feet  with  mean  amplitude  of  5.6  feet.   In 
general,  the  magnitude  of  the  shorter  period  waves  was 
found  to  be  inversely  proportional  to  the  gradients  in 
which  they  were  found.   Theoretically,  internal  waves 
have  periods  ranging  between  the  lower  inertial  and  higher 
Vaisala  or  stability  frequencies.   LaFond  found  periods 
ranging  from  4  to  10  minutes  superimposed  on  longer  diurnal 
cycles  and  four  and  one-half  day  cycles.   He  also  found 
that  internal  waves  moving  towards  the  shore  had  an  average 
speed  of  0.31  knots  and  ranged  from  0.11  to  0.60  knots. 

Physical  intuition  would  lead  us  to  suspect  that  tem- 
perature fluctuations  associated  with  internal  waves  should 
exhibit  a  linear  relationship  between  the  temperature 
gradient  and  the  wave  amplitude.   Therefore  a  reasonable 
expression  for  the  temperature  fluctuations  associated  with 
internal  waves  would  be 


6.   (t) 

1W 


ii 
a  z 


Vt: 
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where  the  overbar  indicates  the  time-averaged  mean  and 

A  (t)  is  the  instantaneous  amplitude  of  the  internal  wave. 
o 

The  temperature  fluctuations  due  to  internal  waves 
are  usually  easily  identified  because  of  the  narrow  range 
of  frequencies  where  significant  energy  is  present,  and 
the  presence  of  a  valley,  in  the  spectral  energy  density, 
separating  the  surface  and  internal  wave  regions.  Further, 
energy  contained  in  the  turbulent  temperature  fluctuations 
is  small  relative  to  the  energy  of  the  internal  waves,  and 
is,  therefore,  only  a  minor  feature  at  these  low  frequen- 
cies. 
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Ill .    EXPERIMENT  AND  DATA  ANALYSIS 

A.   EXPERIMENT 

Measurements  were  made  at  various  depths  with  either 
a  vertical  or  horizontal  orientation  of  the  temperature 
array  as  summarized  in  Table  I. 


TABLE  I.   Summary  of  Experimental  Set-up  for  Minard's 
May  1973  Experiment 


Run 

Number 

Date 

Time 

Depth  to 
top  of 

frame (m) 

Array 

disposition 

Internal 

waves 

Digi- 
tized 

Tide 

1 

15 

May 

1300 

2.7 

vertical 

no 

yes 

ebb 

2 

15 

May 

1437 

5.2 

vertical 

yes(l) 

yes 

slack 

3 

15 

May 

1607 

11.3 

vertical 

yes (2) 

yes 

flood 

4 

15 

May 

1722 

2.7 

vertical 

no 

yes 

flood 

5 

16 

May 

1017 

2.7 

horizontal 

no 

no 

ebb 

6 

16 

May 

1050 

2.7 

horizontal 

yes(l) 

yes 

ebb 

7 

16 

May 

1244 

8.5 

horizontal 

yes (1) 

yes 

ebb 

8a 

16 

May 

1410 

13.1 

horizontal 

no 

yes 

slack 

Figure  5  illustrates  the  array  configuration  for  the 
first  four  runs,  conducted  on  15  May,  where  the  placement 
of  the  thermistors  was  vertical.   Figure  6  illustrates 
the  configuration  of  the  array  on  the  second  day,  16  May, 
when  the  thermistors  were  in  a  horizontal  line.   Run  6 
was  conducted  in  this  configuration  but  in  all  subsequent 
runs  the  array  was  rearranged  such  that  Thermistor  1  was 
removed  and  placed  24  inches  (60.8  cm)  below  Thermistor  5. 
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Figure  5.   Configuration  of  Array  for  Runs  1,  2,  3,  4 
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Figure  6 


Configuration  of  Array  for  Run  6.   Tor  Runs 
7  and  8a,  Thermistor  1  was  removed  and  placed 
60.8  cm  below  Thermistor  5. 
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This  was  done  to  give  an  indication  of  the  vertical  tem- 
perature gradient.   Data  were  recorded  on  a  Sangamo  3562 
analog  tape  recorder  using  F.M.  electronics.   These  data 
were  played  back  at  60  ips  (32  times  speed-up)  and  digi- 
tized at  NPS  on  a  hybrid  (analog-digital)  computer,  CI5000 
/XDS  9300.   The  sampling  rate  of  the  data  analyzed  was 
reduced  to  6.25  Hz  leading  to  a  Nyquist  frequency  of  3.1 
Hz.   This  reduction  was  not  considered  critical,  since 
the  spectra  became  essentially  flat  and  the  signal  lost 
in  noise  before  the  3.1  Hz  frequency  was  reached.   Data 
were  analyzed  on  the  school's  IBM  360-67  computer.   The 
environmental  conditions  consisted  of  light  winds,  small 
amplitude  waves  and  swell,  low  overcast  skies,  and  occa- 
sional internal  waves  during  the  two  day  period.   Further 
details  of  both  the  experimental  set-up  and  the  analysis 
particulars  can  be  found  in  Minard  (1973). 

B.    DATA  PREPARATION 

Initially  all  the  signals  were  graphically  displayed 
on  an  eight  channel  strip  recorder  for  visual  inspection 
of  the  records,  polarity,  error  and  anomaly  checks.   Con- 
firmation was  made  that  runs  5  and  8  were  unusable  for 
digital  purposes.   Run  9  was  designated  Run  8a  and  re- 
placed Run  8,  for  analysis  purposes,  since  they  were  made 
at  the  same  depth  with  the  same  configuration  and  in 
nearly  the  same  time  frame.   The  analysis  was  carried  out 
using  the  nine-track  digital  tapes  previously  prepared 
from  the  analog  records  by  Minard  (1973) . 
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C.   METHOD  OF  ANALYSIS 

The  variances  and  covariances  were  calculated  from 
the  time  series  of  each  run  by  first  removing  any  long 
period  trend  and  any  D.C.  bias  component.   In  addition, 
some  of  the  time  series  were  digitally  high-pass  fil- 
tered to  filter  out  waves  of  period  greater  than  100 
seconds  prior  to  computation.   Spatial  correlations  were 
computed  by  averaging  the  covariance  (C.V.)  between  the 
first  thermistor  and  each  of  the  others  such  that: 

C.V.  =  —  /Th  (t,x)Th.  (t  ,x  +  Ax±)dt 

where  T  is  the  length  of  record  (=  20  min) .   The  C.V. 
was  then  plotted  against  the  appropriate  separation  dis- 
tance along  the  array  using  the  first  thermistor  as  a 
reference.   The  variances  and  covariances  were  normal- 
ized with  respect  to  the  variance  of  the  first  thermis- 
tor in  the  array.   Both  the  normalized  and  the  unnormal- 
ized  spatial  correlation  functions  were  Fourier 
transformed  to  obtain  the  normalized  and  unnormalized 
wavenumber  spectra  for  the  runs. 

Next,  the  gradients  present  in  Runs  1,  4,  7,  and  8a 
were  calculated  since  they  were  to  be  used  in  the  theo- 
retical wave-induced  temperature  spectrum  calculation. 
First,  the  digitized  signals  had  to  be  re-converted  to 
a  true  temperature  since  during  the  recording  and  digi- 
tizing processes  a  number  of  gains  and  biases  had  been 
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introduced.   The  original  temperature  recording  had  been 
a  voltage  conversion  of  the  temperature  variation  around 
a  balance  resistance  in  the  wheatstone  bridge  circuit  of 
each  thermistor.   Therefore,  the  actual  temperature 
would  have  been  of  the  form 

T  =  T  +  AT  , 

where  the  overbar  indicates  the  mean  temperature  and  AT 
is  the  variation  about  the  mean. 

In  San  Diego,  all  signals  were  amplified  before  re- 
cording.  Therefore  the  recorded  signal  would  be  Gl'AT, 
where  Gl  was  the  appropriate  gain.   During  the  digitizing 
process  a  bias  (Bl)  was  applied  prior  to  the  computer 
gain  (G2) ,  which  was  wired  into  the  analog  circuitry.  The 
digitized  signal,  in  terms  of  the  original  fluctuation, 
would  be 


D.S 


G2  (Gl  *  AT  +  Bl)   . 


The  original  temperature  in  terms  of  subsequent  biases 
and  gains  would  be 


T  =  T  + 


(  D.S .  -  B1*G2) 


( 


G2*G1 


)   ' 


The  gradients  were  then  calculated  for  the  selected 
runs  using  the  temperature  difference  between  thermistors 
1  and  5,  which  had  a  constant  spatial  separation  of  24 
inches  (0.6  m),  as  representative  of  the  whole  array. 
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A  separate  program  was  written  to  perform  the  spec- 
tral analysis  of  the  time  series.   The  data  were  read 
into  files  corresponding  to  Thermistor  5  and  the  Baylor 
Wave  Gauge.   The  records  were  then  detrended  and  the 
mean  removed.   Autocorrelation  functions  for  each  file 
and  a  crosscorrelation  function  between  them  were  calcu- 
lated.  The  autocorrelations  were  normalized  with  res- 
pect to  their  variance  and  the  crosscorrelations  with 
respect  to  the  square-root  of  the  product  of  the  indi- 
vidual file's  variances.   From  these  functions,  the  auto- 
and  cross-spectra  were  calculated  by  Fourier  transforming 
the  respective  temporal  correlation  functions,  after 
application  of  a  Parzen  window  filter  to  account  for  the 
finite  length  of  the  time-series. 

From  the  Baylor  Wave  Gauge  spectrum, the  wave-induced 
temperature  spectrum  was  calculated  using  the  spectral 
transfer  function  derived  in  Equation  (5).   In  addition, 
coherence  and  phase  angle  were  calculated  between  the 
temperature  and  waves.   The  theoretical  turbulence  spec- 
trum was  calculated  from  the  difference  between  the 
actual  temperature  spectrum  and  the  wave-induced  tempera- 
ture spectrum  as  given  in  Equation  (9) . 
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IV.        RESULTS 


A.  ANALOG    TEMPERATURE     RECORDS 

Examination  of  the  analog  records  of  temperature,  re- 
produced in  Figures  7  (Run  6)  and  8  (Run  1) ,  show  that  an 
internal  wave  was  present  in  Figure  7  and  not  in  Figure  8. 
The  common  features  of  each  are  the  high  levels  of  tempera- 
ture fluctuation  and  the  general  consistency  of  prominent 
features  throughout  the  temperature  array. 

B.  TEMPERATURE  FIELDS 

Some  long  term  fluctuations  were  noticed  on  all  runs 
and  four  runs  (2,3,6,7)  show  significant  internal  wave 
activity . 

Runs  1  and  4  show  patches  of  increased  turbulence 
which  may  or  may  not  have  been  remnants  of  internal  waves. 
In  any  case,  the  scale  of  the  turbulence  was  such  that 
the  whole  of  the  vertical  array  was  affected.   Figure  9 
is  a  plot  using  the  digitized  data  of  the  first  15  minutes 
of  Run  1  with  a  data  point  plotted  every  4  seconds.   The 
total  vertical  extent  of  the  array  covered  1.625  meters 
and  over  such  a  short  distance  the  temperature  difference 
between  thermistors  was  never  greater  than  0.5°C  and 
usually  less. 

Figure  10,  a  plot  of  the  isotherms  for  Run  1,  was  made 
by  converting  each  thermistor  signal  to  an  actual 
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temperature  reading  in  the  manner  described  in  Section 
III-B.   Subroutine  "CONTUR" ,  from  the  School's  IBM/360 
source  library,  was  used  to  plot  the  isotherms  at  inter- 
vals of  0.3°C  with  a  data  point  every  4  seconds  of  real 
time.   The  temperature  data  from  the  time  series  was 
linearly  interpolated  every  3  inches  across  the  array, 
at  each  time  increment,  as  input  to  the  subroutine.   The 
plot  shows  a  relatively  warm  layer  for  the  first  15  cm, 
underlain  by  a  relatively  cold  layer  in  the  next  15  cm, 
which  in  turn  is  underlain  by  a  relatively  warm  layer 
covering  the  rest  of  the  array. 

Although  no  measurements  were  made  of  salinity  in  the 
detail  required  to  solve  this  question,  the  salinity 
gradient  must  be  slightly  positive  or  at  least  isohaline 
for  the  first  30  cm  of  the  array,  and  then  increasingly 
positive  to  at  least  60  cm  in  order  for  this  system  to  be 
dynamically  stable.   If  this  is  the  case,  then  the  possi- 
bility of  double— diffusion  type  microstructur e ,  with 
stable  stratification  and  convective  mixing  within  each 
layer,  exists  between  the  central  and  lower  layer.   The 
increased  gradient  of  temperature  between  15  and  30  cm 
is  consistent  with  an  inherently  stable  double-diffusion 
mechanism.   If,  in  fact,  the  salinity  gradient  is  negative 
in  the  first  30  cm  then  the  possibility  of  salt-fingering 
exists  (inherently  unstable  double-diffusion),  although 
there  is  no  direct  evidence  in  the  isotherm  field  to  sug- 
gest this  mechanism  is  occurring. 
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The  decrease  in  the  temperature  gradient  in  the  lower 
half  of  the  array  (Figure  11)  may  be  an  artifact  produced 
by  the  array  spacing  and  the  subsequent  interpolation  of 
points  in  the  contour  plot. 

The  blobular  formations  within  the  warm  upper  layer 
and  in  the  cold  central  layer  are  probably  due  to  inter- 
actions with  small  amplitude  internal  waves  which  were  not 
recognized  as  such  from  the  temperature  records  alone.   It 
would  appear  that  the  layers  are  being  alternately 
squeezed  and  relaxed  with  a  periodicity  approximating 
3  minutes.   These  blobs  compare  well  in  time  to  the  in- 
creases and  decreases  in  the  turbulent  temperature  evident 
in  Figure  9. 

Runs  which  show  significant  internal  wave  activity  are 
itemized  in  Table  II  by  run  number  and  "periodicity" 
measured  from  the  analog  records.   Periodicity  is  at  best 

TABLE  II.   Internal  Wave-Periodicity 


run  number 

period  ( sec ) 

2 
3 
6 
7 

1350 
1050 
1450 
1300 

approximate  due  to  the  limited  number  of  internal  waves 
per  record.   Figure  11  is  a  plot  of  the  first  15  minutes 
of  Run  3  showing  the  temperature  variation  for  Thermistor  1 
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(top  of  the  array,  at  11.3  meters  below  the  surface)  and 
Thermistor  6  (towards  the  bottom  of  the  array) .   One  in- 
ternal wave  had  just  finished  at  the  start  of  the  plot. 
A    second  was  starting  at  the  end  of  the  plot  and  con- 
tinues into  Figure  12,  the  next  15  minutes  of  Run  3.   It 
illustrates  the  periodicity  of  the  internal  wave.   Figure 
13  is  the  computer  plot  of  the  isotherms  for  the  first  15 
minutes  of  Run  3.   The  intrusion  of  colder  water  into  the 
upper  portions  of  the  array  and  its  subsequent  withdrawal 
as  the  internal  wave  passes  is  quite  noticeable.   The 
salinity  profile  corresponding  is  unavailable  but  is 
probably  a  positive  or  isohaline  gradient  with  increasing 
depth.   Double-diffusion  type  micros tructure  would  not  be 
present  under  these  conditions  of  stratification  and 
billow  turbulence  microstructure  would  be  more  likely. 
The  thermistor  array  during  Run  3  was  located  just 
below  a  slight  thermocline  which  had  been  decreasing  in 
depth  and,  to  some  extent,  intensity  from  the  time  of 
Run  1 (Minard  1973-Figure  1) .   This  suggests  that  the 
presence  of  internal  waves  in  the  period  of  time  from 
Run  1  to  Run  3  was  contributing  to  the  depression  and 
erosion  of  the  thermocline.   Figure  13  shows  the  extreme 
variability  and  fluctuation  in  an  individual  isotherm. 
Some  of  the  fluctuation  is  no  doubt  attributable  to  tur- 
bulence which  could  be  a  cause  of  the  thermocline  erosion 
However,  the  assumption  that  the  fluctuations  are  solely 
caused  by  turbulence  must  be  tempered  to  some  degree 
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since  successive  peaks  in  the  fluctuations  are  separated 
by  about  12  seconds,  a  figure  that  is  suspiciously  close 
to  the  surface  wave  period. 

C.    SPATIAL  CORRELATIONS 

Table  III  contains  a  listing  of  the  covariances  be- 
tween Thermistor  1  and  the  others,  and  the  individual 
variances  of  each  thermistor,  for  all  of  the  runs  for  both 
non-filtered  and  high-pass  filtered  data.   To  accomplish 
the  high-pass  filtering,  a  Fourier  analysis  was  made  of 
each  20  minute  segment  of  data;  the  Fourier  coefficients 
generated  were  used  to  generate  a  time  series  describing 
the  waves  of  100  seconds  and  greater.   This  series  was 
then  subtracted  from  the  original  time  series  which  was 
then  the  filtered  time  series  used  for  analysis. 

One  of  the  most  common  ways  to  define  a  scale  length 
of  some  inhomogeneity  is  to  measure  the  length  required 
for  the  correlation  to  fall  to  1/e  (0.368)  of  its  original 
value.   The  correlation  functions,  shown  in  Figures  14  and 
15,  are  of  a  very  long  scale  feature  and  only  one  run  in 
each  figure  actually  falls  to  0.368.   This  phenomena  is 
interpreted  as  being  due  to  the  extremely  long  scale  of 
the  internal  waves  and  their  attendant  higher  overall 
energy.   Consequently  it  v/as  decided  to  filter  out  the 
waves  of  periods  100  seconds  and  longer.   This  was  accom- 
plished digitally,  as  explained  earlier.   The  results, 
after  the  digital  filtering  of  the  temperature  signal, 
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are  illustrated  in  Figures  16  and  17  where  it  is  appar- 
ent that  waves  of  100  sec  and  longer  had  been  contribut- 
ing substantial  energy  as  evidenced  by  the  overall 
lowering  of  the  spatial  correlation  curves  when  filtered, 

Table  IV  lists  the  correlation  lengths  for  each  run 
after  the  original  signal  had  been  high-pass  filtered. 

TABLE  IV.   Correlation  Lengths 


Run 

Depth-Surface 

to 

Correlation 

Number 

Top  of  Array 

(m) 

Length  (cm) 

1 

2.7 

100.  5 

2 

5.2 

100.  5 

3 

11.  3 

35.  2 

4 

2.7 

130.  5 

6 

2.7 

40.0 

7 

8.5 

50.  0 

8a 

13.  1 

21.4 

Although  it  might  be  expected  that  the  shortest  cor- 
relation lengths  would  occur  at  the  surface,  where  waves 
dominate  the  structure,  no  such  relationship  is  apparent. 
Perhaps  this  can  be  explained  by  the  fact  that  the  array 
dimensions  are  small  compared  to  the  depth  and  breadth  of 
the  surrounding  ocean  area;  or  the  scale  of  the  turbulence 
may  be  quite  large  compared  to  the  depth  of  the  water;  or 
maybe  it  further  illustrates  that  waves  do  little  mixing 
but  only  pump  the  water  up  and  down. 

A  consistent  trend  was  noticed  in  that  the  horizontal 
array  always  had  shorter  correlation  lengths  than  the 


53 


in 
o 


CN 


0 

W 

2      •      » 
<  ^    in 

dm   n 


H     •      • 
Pl,  cn    in 
w 
a 


z    z 

a:    oc 

1 

< 
I 


ro 
n 
m 


Z 

I 

o 

I 


e 
o 


0) 

o 
c 

r0 
■P 
CO 

■rH 

Q 

0> 

(0 


>i 
< 


ra 
u 

•H 

-P 
U 
0) 

> 


c 
o 

■H 
-P 
ffl 
cH 
Q) 
U 
U 
O 

u 


rd 

•H 
4-> 
rd 
& 


<D 
U 

P 

•H 


«-    UOT^PXOJ[J03    pOZTTPUUON 


54 


in 


o 

•-t 

X 

o 

o 

*■* 

w 

CJ 

• 

• 

CM 

^ 

ri* 

en 

a\ 

H 

n 

ro 

05 

VD 

^3* 

Cn     eg     CD 


Z  Z 

D  D 

OH  DC 

O  ' 
I 


o 

CO 

Z 

« 
I 

a 
i 


0) 

o 

c 
<d 
-p 
w 

•H 
Q 

en 

(t! 

*1 


>1 


■P 
o 

N 
•H 

5-1 
O 


o 
•p 

iH 

(1) 
H 
M 
O 
U 


■P 


0) 
•H 


UOT^T?T3JJ03    pDZTTICUIJON 


vertical  array  at  similar  depths.   However,  this  result 
is  not  conclusive  since  horizontal  and  vertical  array 
runs  were  made  on  different  days. 

Further  examination  of  Table  III  shows  that  the  runs 
with  significant  internal  wave  activity  have  higher 
energy  as  evidenced  by  the  higher  numerical  values  of 
variance  and  covariance.   This  trend  is  present  even 
when  the  signal  was  filtered  for  waves  of  100  sec  and 
greater  except  that  Runs  2  and  3  are  not  significantly 
greater  in  energy  than  Runs  1,  4,  and  8a.   This  indicates 
that  Runs  6  and  7  contain  significantly  more  energy  at 
all  frequencies  than  Runs  2  and  3.   Also  the  difference 
in  covariances  between  signals  not  filtered  and  filtered 
show  that  Runs  1  and  8a  contained  much  less  energy  at 
periods  greater  than  100  sec  than  any  of  the  other  runs. 

D.   WAVENUMBER  SPECTRA 

The  normalized  wavenumber  spectra  were  calculated  by 
Fourier  transforming  the  spatial  correlation  functions 
(Figures  16  and  17)  calculated  as  explained  in  Section 
III.B.   The  spatial  correlation  was  calculated  at  equal 
incremenr.s  along  the  length  of  the  array  by  linearly  in- 
terpolating between  the  calculated  points  on  the  array. 
Figure  18  is  a  composite  of  the  normalized  wavenumber 
spectra  when  the  array  v/as  vertical  (15  May)  ;  and  Figure 
19  is  a  composite  when  the  array  was  horizontal  (16  May). 
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Figure  18.   Vertical  Array  Wavenumber  Spectra 
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The  spectra  show  a  high  degree  of  consistency  between 
runs.   In  the  vertical  array,  Runs  2  and  3  show  slightly 
more  energy  at  all  wavenumbers  than  Run  4  but  less  energy 
than  Run  1.   It  would  appear  that  Run  4  had  a  major  por- 
tion of  its  energy  at  periods  greater  than  100  sec  as  was 
the  case  for  Runs  2  and  3  with  the  internal  waves.   In 
the  horizontal  array  spectra,  all  runs  are  very  similar 
in  shape  whereas  the  vertical  array  spectra  varied  con- 
siderably in  shape.   It  is  felt  that  this  may  be  inter- 
preted as  a  difference  in  the  source  of  turbulence  on  the 
two  days . 

Each  of  the  spectra  has  what  is  believed  to  be  an 
anomalous  spectral  energy  increase  at  higher  wavenumbers. 
Gregg,  Cox  and  Hacker  (1973)  observed  a  similar  increase 
at  slightly  higher  wavenumbers  in  the  temperature  gradi- 
ent spectra.   They  attributed  the  initial  rapid  fall-off, 
below  0.06  radians  per  meter  (rpm),  to  dynamical  processes 
dominated  by  the  nearly  exponential  mean  temperature 
gradient  with  depth.   For  wavenumbers  greater  than  0.06 
but  less  than  31.4  rpm  they  felt  the  spectrum  was  domi- 
nated by  the  irregularly  spaced  high  gradient  regions.  In 
the  region  of  interest,  above  31.4  rpm,  where  high  levels 
of  micros tructure  exist  the  spectrum  takes  the  form  of  a 
minimum  followed  by  a  maximum  at  higher  wavenumbers.   This 
last  region  is  in  general  the  same  wavenumber  space  as  the 
increase  in  our  spectrum.   Whether  the  increase  is  due  to 
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the  same  cause  is  uncertain;  but  there  were  high  levels 
of  microstructure  activity  in  all  the  runs.   However,  it 
may  be  that  this  increase  is  due  to  other  factors  such 
as  the  manner  in  which  the  correlation  function  was  calcu- 
lated in  deriving  the  spectrum.   No  attempt  was  made  to 
fit  the  correlation  function  to  other  interpolation  forms, 
e.g.,  Gaussian  or  exponential,  near  the  origin  since  the 
structure  function  when  calculated  did  not  show  that  the 
time  series  were  particularly  stationary. 

It  seems  reasonable  that  the  wave-induced  contribu- 
tion to  the  wavenumber  spectra  would  occur  over  a  narrow 
frequency  band  corresponding  to  the  frequencies  of  high 
coherence  in  the  frequency  spectra.   Consequently,  the 
spatial  scale  in  wavenumber  space  should  show  corres- 
pondence to  the  orbital  excursion  of  the  wave-induced 
particle  motion  within  this  narrow  band.   The  scale  of 
the  particle  displacement  is  calculated  using  the  root- 
mean-s^quare  (rms)  wave  height  in  each  run.   The  vertical 
displacement,  £(t),  using  linear  wave  theory,  is  related 
to  the  wave  profile,  n.(t),  by 


=  [s1nh_Mz1d)_] 
si nh  kd 


where  the  term  in  brackets  is  the  spectral  response  func-' 
tion,  H  (a).  The  maximum  displacement  can  be  approximated 
at  the  peak  in  the  frequency  spectrum  as 


max 


H  (a)    n 

w    max  rms 

ae/8z 


(12) 


6  0 


where  H  (a)     is  the  response  function  at  the  frequency 
w    max 

of  maximum  coherence  and  n     is  the  rms  wave  height. 

rms 

The  horizontal  excursions  can  be  similarly  calculated, 
using  linear  wave  theory,  as 


f (t)  =  [£    ]  coth  k(z+d) . 
max 


The  results  of  these  calculations  are  presented  in  Table 
V. 


TABLE  V.   Particle  Displacements 


Run 

Wave 

Height 

(m) 

Depth 
of  Th. 
5  (m) 

Freq. 
(Hz) 

max 
(m) 

Wave- 
num- 
ber 
(rad/m) 

f 

max 

(m) 

Wave- 
num- 
ber 
(rad/m) 

Expected 
wave- 
number 
(rad/m) 

1 

0.99 

3.6 

0.092 

0.76 

8.3 

1.25 

5.0 

9.0 

6 

1.27 

3.6 

0.092 
0.214 

0.98 
0.66 

6.4 
9-5 

1.63 
0.66 

3.9 

9.5 

8.0 
8.0 

7 

1.52 

9.3 

0.092 

0.68 

9.2 

1.63 

3.9 

6.0 

8a 

; 

1.38 

13.8 

0.092 

0.30 

21.1 

1.41 

4.5 

18.0 

These  scales  are  representative  of  the  diameter  of  the 
water  particle  orbit.   It  would  be  expected  that  if  there 
is  a  peak  in  the  frequency  spectrum  (Figure  24)  corres- 
ponding to  surface  wave-induced  motion,  then  there  should 
be  a  peak  in  the  wavenumber  spectrum  associated  with  the 
same  motion  at  the  wavenumber  of  the  excursion.   Examina- 
tion of  Figures  18  and  19  shows  that  although  there  is  no 
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peak  in  the  spectra  there  is  a  definite  change  in  the 
slope  of  each  spectrum  in  the  vicinity  of  the  scale 
length  calculated  through  Equation  (12).   The  fact  that 
the  intersections  do  not  lie  exactly  at  the  change  of 
slope  may  be  due  to  the  use  of  the  rms  wave  height 
rather  than  the  median  height  or  some  other  value.  Also 
as  would  be  expected,  the  scale  length,  E, ,    decreases  as 
the  depth  increases  illustrating  the  decrease  of  the 
orbital  motion  with  depth. 

E.    GRADIENT  FIELDS 

Gradients  for  Runs  1,  4,  7,  8a  were  calculated  and 
the  digital  data  plotted.   Table  VI  is  a  compilation  of 
the  calculated  mean  gradient,  the  coherent  gradient  used 
in  calculating  the  wave-induced  spectrum,  and  the  gradi- 
ents as  determined  by  Minard  from  the  bathythermographs 
[Minard,  1973]  . 

TABLE  VI.   Gradient  Data 


Top  of  the 

Top  of  the 

B.T. 

Calculated 

Coherent 

Run 

thermo- 

frame 

gradient 

gradient 

gradient 

Number 

cline  (m) 

(m) 

(°C/m) 

(°C/m) 

(°C/m) 

1 

4.25 

2.7 

0.98 

0.93 

.009 

4 

14.7 

2.7 

0.25 

0.19 

not  calc. 

6 

0 

2.7 

0.44 

not  calc. 

.096 

7 

14.5 

8.5 

0.11 

-0.41 

-.07 

8a 

5.0 

13.1 

0.54 

0.10 

.017 

G2 


As  can  be  seen  from  Table  VI,  the  gradients  taken 
from  the  Bathythermograph  (BT)  are  in  general  agreement 
with  the  time-averaged  mean  gradient  from  the  time-series, 
with  the  exception  of  Run  7.   The  difference  can  be 
accounted  for  because  of  the  difference  in  length  over 
which  the  gradients  were  determined.   The  gradient  deter- 
mined from  the  BT's  was  originally  averaged  over  a  dis- 
tance of  18  meters  and  was  all  that  was  available,  whereas 
the  calculated  mean  gradient  was  determined  over  24  inches 
(0.6  m) ,  the  distance  between  Thermistors  1  and  5  on  all 
runs.   The  difference  in  length  over  which  the  gradient  is 
calculated  is  important  and  can  be  illustrated  from  the 
results  of  Run  1.   The  mean  calculated  gradient  was  deter- 
mined to  be  0.93°C/m  (Table  VI).   Figure  9  shows  the  tem- 
perature profiles  of  Thermistors  1  and  6  of  Run  1  (same 
time  frame  but  Thermistors  1  and  6  are  separated  by  14 
inches  (0.30  m)  more  than  Thermistors  1  and  5) .   Here  the 
maximum  temperature  difference  is  about  0.5°C,  and  since 
the  separation  was  38  inches  (0.96  m)  between  thermistors, 
the  gradient  would  be  a  maximum  of  0.52°C/m.   It  can  be 
seen  from  this  that  a  small  change  in  the  length  can  re- 
sult in  a  large  change  in  the  gradient. 

It  can  also  be  seen  that  the  BT  and  calculated  gradi- 
ents are  much  greater  than  the  gradient  required  to  match 
coherences  in  the  wave-induced  spectrum.   This  indicates 
that  the  spectral  model  would  overestimate  the  wave-induced 
temperature  fluctuations  if  the  actual  mean  gradient  were 
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used.   Thornton  et  al  (1974)  noticed  the  same  result 
and  calculated  the  effect  of  buoyancy  to  try  and  account 
for  some  of  the  difference.   They  concluded  that  "buoy- 
ancy effects  do  not  dampen  the  motion,  but  serve  to 
amplify  it  relative  to  the  homogeneous  case."   Hence 
they  could  not  explain  why  the  mean  gradients  are  greater 
than  the  coherent  gradients  since  the  wave  solution  would 
be  improved  by  less  than  1%  by  including  the  buoyancy 
effects.   What  is  obvious  is  that  the  specification  of 
a  constant  mean  gradient  is  not  completely  realistic  in 
the  presence  of  internal  waves  and  turbulence. 

Figure  20  is  a  plot  of  the  digitized  data  of  Thermis- 
tors 1  and  5  during  the  first  22  minutes  of  Run  1. 
Thermistor  1  was  0.6  m  above  Thermistor  5  and  therefore 
a  positive  gradient  is  indicated,  i.e.,  temperature  de- 
crease with  depth.   Figure  21  shows  both  the  gradient  and 
Thermistor  number  1  temperature  plotted  together  for  the 
first  22  minutes  of  Run  1.   Each  point  on  the  plot  is  the 
averaged  gradient  over  1.44  seconds.   Thermistor  1  is 
replotted  from  Figure  20  but  at  a  slightly  different 
scale.   The  arrows  labelled  A  and  B  in  both  Figures  20  and 
21  serve  to  indicate  noticeable  points  of  gradient  increase 

Figures  22  and  23  show  the  same  information,  as  ex- 
plained above,  for  Run  7  where  Thermistor  5  is  now  0.6  m 
above  Thermistor  1.   Again  a  point  is  marked  (C)  to  illus- 
trate an  increase  in  gradient.   What  is  interesting,  how- 
ever, is  the  difference  in  the  temperature  profiles  and 
gradients  between  Runs  1  and  7. 
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Figure  20.   Digital  Temperature  Record,  Run  1 
Thermistors  1  and  5. 
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Figure  21.   Temperature  Gradient  Plot,  Run  1.   Gradient 
between  Thermistors  1  and  5,  and  Temperature 
Profile  of  Thermistor  1. 
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In  Run  1,  the  temperature  profile  appears  to  be  as 
'turbulent'  and  perhaps  more  so  than  Run  7.   Run  7,  how- 
ever, is  marked  at  the  beginning  by  the  passage  of  an 
internal  wave.   The  gradients  in  Run  1  are  relatively 
constant,  being  uniformly  positive  with  a  deviation  of 
less  than  0.5°C/m  from  the  mean.   On  the  other  hand,  the 
gradients  in  Run  7  are  extremely  variable ,  being  both 
positive  (indicating  stable  conditions)  and  negative 
(indicating  unstable  conditions)  with  extremes  ranging 
from  +0.58  to  -1.04°C/m.   In  addition,  it  appears  that 
most  of  this  record  is  best  characterized  as  dynamically 
unstable.   This  is  certainly  indicative  of  intense  turbu- 
lence associated  with  the  internal  wave  and  there  may  be 
some  indication  that  the  wave  has  'broken'.   In  the  area 
of  point  D  on  Figure  23,  the  gradient  is  almost  zero  for 
a  period  of  nearly  150  sec  and  represents  a  zone  of 
nearly  homogeneous  water  brought  about  from  the  intense 
turbulent  mixing.   In  an  analog  to  surface  waves,  this 
zone  could  be  termed  an  'internal  swash  zone'.   The  calcu- 
lated gradients  for  Runs  4  and  8a  are  similar  to  the  re- 
sults for  Run  1. 

F.    TEMPERATURE  FREQUENCY  SPECTRA  ANALYSIS 

The  results  of  Run  6  will  be  discussed  more  fully 
than  the  other  runs  since  it  most  clearly  demonstrates 
the  effect  of  surface  wave-induced  fluctuations.   The 
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other  three  runs  analyzed  (1,  7,  8a)  will  be  discussed 
only  insofar  as  they  differ  from  Run  6. 

1 .   Surface  Wave-Induced  Temperature  Fluctuations 

Equation  (7)  was  used  to  calculate  the  wave-in- 
duced temperature  spectra.   Initially  the  mean  tempera- 
ture gradient,  as  calculated  from  the  time-series,  was 
used  in  the  transfer  function  given  in  Equation  (5). 
This  led  to  an   overestimation  in  the  wave-induced  tem- 
perature spectra.   The  variability  in  the  temperature 
gradient  field  was  considered  to  be  mainly  a  function  of 
the  internal  waves  present  and  consequently  the  use  of  a 
mean  gradient  was  hardly  realistic.   If  an  approximation 
to  a  mean  gradient  had  to  be  used,  it  may  have  been 
better  to  use  the  median,  or  most  commonly  occurring, 
gradient.   Since  a  statistical  analysis  of  the  gradient 
was  not  possible  it  was  decided  to  use  the  method  of 
establishing  a  coherent  gradient  as  Thornton  et  al  (1974) 
had  done.   Therefore,  a  gradient  was  used  which  matched 
the  measured  peak  coherence  value  between  waves  and  tem- 
perature with  a  calculated  (theoretical)  coherence  using 
the  calculated  wave-induced  temperature  spectrum  substi- 
tuted into  Equation  (11) .   Their  rational  for  doing  this 
was  that  if  the  turbulent  and  wave-induced  temperature 
fluctuations  are  statistically  independent  then  the 
measured  coherence  should  be  closely  approximated  by  this 
equation  provided  the  transfer  function  was  correct. 
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Since  the  model  is  based  on  linear  theory  the  only  "free 
variable"  in  the  transfer  function  is  the  gradient.  This 
coherent  gradient  (Table  VI)  is  significantly  different 
from  the  mean  gradient.   Table  VII  shows  the  appropriate 
values  of  depth,  frequency,  coherence,  temperature  spec- 
tral density,  wave-induced  spectral  density,  and  average 
phase  angle  (between  the  wave  and  temperature  spectra) 
at  the  point  of  maximum  coherence  in  the  wave  and  tempera- 
ture spectra. 

TABLE  VII.   Frequency  Spectra  Data 


Temp. 

W-I  Temp. 

Average 

Depth 

Freq. 

Coherence 

Spectral 

Spectral 

Phase 

Run  # 

(m) 

(Hz) 

ma'x. 

Value 

Value 

Difference 

1 

3.6 

.092 

.23 

,18xl0~4 

,41xl0~5 

145 

6 

3.6 

.092 

.42 

,14xl0~2 

.59xl0~3 

155 

.214 

.44 

.39xl0~3 

.11x10" 

155 

7 

9.3 

.092 

.40 

.62xl0~3 

.25xl0~3 

140 

8a 

13.8 

.09*2 

.37 

.64xl0~5 

.24xl0~5 

105 

Figure  24  shows  the  measured  temperature  and  calculated 
wave-induced  temperature  spectra,  and  the  coherence  and 
phase  angle  between  waves  and  temperature  which  are  dis- 
cussed in  the  following  sections. 
2 .   Coherence 

The  measured  coherence  between  waves  and  tempera- 
ture is  shown  in  the  middle  graph  of  Figure  24  as  the 
solid  line.   The  (theoretical)  coherence  between  waves 
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perature Spectra,  Run  6. 
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and  wave-induced  temperature  was  calculated  from  Equation 
(11)  assuming  statistical  independence  and  is  illustrated 
as  the  dashed  line  on  the  same  figure.   It  is  only  in  the 
relatively  narrow  band  from  0.05  to  0.3  Hz  that  there  is 
significant  coherence  between  the  spectra  of  waves  and 
temperature.   This  corresponds  to  periods  of  between  3.3 
and  20  seconds  which  is  most  certainly  characteristic  of 
surface  waves  and  swell.   The  theoretical  coherence, 
which  was  adjusted  to  match  the  measured  coherence  at  the 
first  peak  value,  approximates  the  measured  coherence 
very  closely.   Hence  the  separation  of  the  two  spectra 
can  be  interpreted,  from  the  increasing  lack  of  coherence 
with  increasing  frequency,  as  an  increasingly  high  ratio 
of  turbulence   (noise)  to  coherent  surface  wave-induced 
temperature  fluctuations. 
3  .   Phase  Difference 

The  last  graph  in  Figure  24  is  the  measured  spec- 
tral phase  difference  between  the  waves  and  the  tempera- 
ture.  A  visual  inspection  of  the  figure  reveals  that 
over  the  range  of  significant  coherency,  the  average  phase 
difference  is  155°  (values  for  the  other  runs  are  listed 
in  Table  VII) .   In  theory,  the  phase  of  the  waves  should 
lead  the  temperature  by  180°.   The  difference  was  un- 
expected since  the  results  of  Thornton  et  al  (1974)  indi- 
cated a  much  closer  relationship  between  theory  and  actual 
measurement.   In  addition  they  found  that  the  span  over 
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which  the  phase  difference  was  close  to  180°  was  much 
wider  than  the  frequency  band  of  high  coherence.   This 
occurred  in  these  results  only  once,  in  Run  7.   Also, 
the  phase  difference  appeared  to  be  depth  dependent  be- 
cause in  Runs  6  to  8a  where  the  depth  of  the  sensor 

(Thermistor  5)  increased,  the  phase  difference  decreased 
from  155  to  105°. 

The  reasons  for  these  differences  are  not  readily 
forthcoming.   In  Run  1,  as  noted  earlier,  the  stratifica- 
tion was  different  from  the  other  runs  in  that  over  the 
length  of  the  array  the  temperature  field  was  highly  sug- 
gestive of  some  double-diffusion  mechanism  operating. 
This  may  suggest  that  in  such  cases  buoyancy  plays  a 
much  greater  part  in  the  destruction  of  coherence  between 
waves  and  temperature  than  in  cases  where  temperature 
steadily  decreases  with  depth. 

The  actual  difference  between  the  average  phase 
difference  and  180°  is  even  more  puzzling.   With  their 
sensor  at  5.5  m  Thornton  et  al  (1974)  found  almost  exact- 
ly a  180°  phase  difference  from  about  .05  to  0.22  Hz. 
The  results  of  this  experiment  at  a  lesser  depth  indicate 
a  phase  change  of  about  155°  (about  25°  difference)  which 
decreased  as  the  sensor  depth  increased.   At  first  the 
difference  was  attributed  to  the  horizontal  displacement 

(leading  to  a  phase  lag)  of  Thermistor  5  from  the  Baylor 
Wave  Gauge.   The  first  order  expressions  for  the  wave  and 
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temperature  fluctuations  are: 


l"|  =  a  cos  (kx  -  at ) 


6  =  H  (9)  a  cos  (kx  -  at) 
a 


■H  (0)  a  cos  [180°  -  (kx  -  at)]  . 


Using  the  wave  gauge  as  a  reference  position,  the  phase 
angles  become, 

I   =  -at 

n 

£„  =  180°  -  (KAx  -  at)  , 
y 

so  that  the  phase  difference,  AE  =  EA  -  £  ,  is  180°  if  Ax 
is  zero.   The  horizontal  displacement  was  0.19  m  which 
for  this  calculation  at  the  frequency  of  maximum  coherence 
yielded  a  difference  of  0.52°  from  180°.   This  is  a 
definite  indication  that  all  is  not  as  simple  as  the 
assumptions  in  the  model  would  have  us  believe.   Internal 
waves,  which  of  themselves  would  be  180°  out  of  phase  with 
the  temperature  fluctuations,  would  certainly  complicate 
the  picture  because  of  their  interaction  with  the  surface 
waves  and  the  vertical  displacement  of  the  thermistor  and 
the  wave  gauge.   Buoyancy  effects  brought  about  by  small 
scale  instabilities  such  as  those  noted  in  Run  7  may  be 
another  factor.   Whatever  the  cause,  it  can  be  said  from 
these  results  that  the  phase  difference  is  not  due  to  the 
surface  wave-induced  motions  alone. 
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4 .   Turbulent  Temperature  Fluctuations 

If  in  the  consideration  of  the  theory  we  assume 
that  turbulent  and  wave-induced  temperature  fluctuations 
are  statistically  independent,  then  the  area  between  the 
spectral  energy  density  curves  on  Figure  24  should  repre- 
sent the  turbulent  energy  as  given  by  Equation  (10) . 
Figure  25  is  a  plot  of  the  turbulent  temperature  spectrum 
for  Run  6.   Thornton  et  al  (1974)  calculated  the  turbulent 
temperature  spectra  considering  the  possibilities  that 
turbulent  and  wave-induced  temperature  fluctuations  were 
correlated  (Equation  (8))  as  well  as  statistically  inde- 
pendent (Equation  (10)).   They  reported  almost  identical 
results.   Hence,  only  the  latter  is  shown  here. 

The  turbulent  spectra  for  all  the  runs  were 
generally  consistent  and  compared  well  with  the. -5/3  slope 
illustrated  in  Figure  25.   There  is  a  definite  indication 
of  a  change  over  in  the  inertial  subrange  at  a  period  of 
about  9  seconds  (f  =  .11  Hz) .   This  is  to  be  expected 
since  turbulence  occurs  at  scales  smaller  than  those 
associated  with  the  waves  from  which  turbulent  energy  is 
derived.   In  fact,  the  slope  of  the  spectra  is  much  nearer 
to  -2  than  -5/3„   Black  (1965)  has  suggested  this  slope  is 
to  be  expected  more  often  than  -5/3  when  there  is  density 
stratification.   This  would  agree  with  the  results  sug- 
gested by  the  calculated  structure  function.  These  results 
appear  to  generally  validate  the  assumption  of  statistical 
independence  between  turbulence  and  wave-induced 
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temperature  fluctuations  and  that  the  natural  preference 
in  nature  is  for  a  near  -5/3  dependence  of  turbulent 
spectral  energy  on  frequency  in  the  inertial  subrange. 

Seitz  (1971)  has  suggested  that  in  a  stationary, 
homogeneous,  isotropic  field  of  turbulence  the  energy 
spectra  will  have  a  dependence  on  a  -2  power  law  at  low 
frequencies  and. a  -5/3  power  law  at  high  frequencies  pro- 
viding convection  velocities  are  moderate.   The  transi- 
tion from  one  power  law  to  the  other  should  be  at  some 
intermediate  frequency  and  would  be  a  function  of  the 
magnitude  of  the  convection  velocity.   Following  this 
line  of  reasoning  then  it  can  be  seen  in  Figure  25  that 
all  runs  exhibit  inertial  subrange  character  in  the  energy 
spectra  above  0.1  Hz„   Runs  1  and  8a,  which  had  no 
apparent  internal  wave  activity,  exhibit  similar  slopes 
of  -7/3  below  0.031  Hz  with  transition  between  these  fre- 
quencies.  However,  Runs  6  and  7,  which  did  show  signifi- 
cant internal  wave  activity,  exhibit  slopes  of  -3  and  -4.5 
respectively  below  0.031  Hz.   This  would  tend  to  reinforce 
the  notion  that  internal  waves  contain  much  more  energy 
at  the  lower  frequencies  and  are  primarily  responsible  for 
the  time  series  being  non-stationary. 
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V.   CONCLUSIONS 

The  temperature  field  is  extremely  complex  in  the 
shallow  waters  near  the  NUC  Tower.   There  were  indica- 
tions that  double-diffusion  type  microstructure  was 
present  as  well  as  the  more  normal,  for  this  area,  billow 
turbulence  microstructure  which  appeared  to  be  suggested 
by  the  temperature  stratification  for  the  majority  of  the 
runs. 

Temperature  fields  with  no  apparent  internal  waves 
may  contain  turbulent  oscillations  suggestive  of  internal 
waves  or  their  remnants.   The  basis  for  this  statement  is 
that  the  oscillations  in  Run  1  appeared  to  be  periodic  of 
period  three  minutes  which  is  close  to  that  reported  by 
LaFond  (1966)  for  internal  waves. 

Significant  internal  waves  present  on  the  two  days  of 
the  experiment  were  periodic  with  approximately  20  minute 
periods.   In  form  they  were  more  like  solitary  waves  rather 
than  sinusoids  but  this  is  probably  due  to  the  fact  that 
they  were  shallow  water  waves  at  this  location.   Judging 
from  the  temperature  profiles  the  internal  waves  appeared 
to  have  positive  displacement  from  a  mean  position  with 
little  or  no  negative  displacement.   The  positive  displace- 
ments occupied  less  than  one-half  of  one  period. 

Internal  waves  appeared  to  be  the  major  source  of  tur- 
bulence during  this  experiment  under  the  conditions  of 
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light  winds  and  small  ampitude  waves.   This  conclusion 

is  based  on  the  higher  variances  in  the  runs  with  internal 

waves . 

When  internal  waves  were  present  there  was  signifi- 
cantly more  energy  at  all  wavenumbers  and  not  only  at  low 
wavenumbers.   This  statement  is  predicated  on  the  fact 
that  even  after  filtering  out  waves  of  period  100  seconds 
and  greater  the  runs  with  internal  waves  still  had  more 
energy  than  runs  without  internal  waves. 

Internal  waves  may  possibly  break,  thus  contributing 
a  large  amount  of  turbulent  energy  to  a  localized  area 
and  thus  assist  mixing. 

With  respect  to  spatial  correlation  lengths,  the  length 
of  the  array  (64  inches  =  1.6  meters)  was  too  short  to 
measure  the  scale  length  of  the  internal  waves.   The  high- 
pass  filtered  time  series  had  correlation  lengths  of  130  cm 
and  less.   There  did  not  appear  to  be  any  depth  dependence 
on  the  length  but  there  may  have  been  a  dependence  on 
whether  the  array  was  vertical  or  horizontal. 

With  respect  to  the  gradient  field,  it  is  confirmed 
that  temperature  gradients  taken  over  short  distances  are 
extremely  variable  over  time  and  space  as  expected  for  a 
turbulent  phenomenon. 

Wavenumber  spectra  may  be  significantly  altered  at 
higher  wavenumbers  depending  on  the  type  of  interpolation 
used  in  determining  the  spatial  correlation  function  to  be 
Fourier  transformed.   Further,  the  vertical  and  perhaps 
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the  horizontal  particle  displacements  at  the  frequency 
of  maximum  coherence  in  the  frequency  spectra  appear  as 
changes  of  slope  in  the  wavenumber  spectra. 

The  conclusions  reached  in  respect  -to  the  spectral 
model  proposed  by  Thornton  et  al  (1974)  were  that: 

1.  The  model  is  not  always  realistic  in  the  speci- 
fication of  a  constant  mean  gradient.   The 
gradient  required  to  be  used  in  the  transfer 
function  to  achieve  consistency  between  the  actual 
temperature  spectrum  and  the  wave-induced  tempera- 
ture spectrum  was  much  less  than  the  observed  or 
calculated  mean  gradients. 

2.  The  wave-induced  temperature  spectrum  of  the 
model  showed  great  similarity  to  the  actual  tem- 
perature spectrum  over  a  narrow  band  of  frequen- 
cies associated  with  regions  of  high  coherence 
between  surface  waves  and  temperature  spectra. 

3.  The  model  describes  reasonably  well  the  wave- 
induced  temperature  fluctuations  in  shallow  water 
(less  than  18  m) . 

4.  The  model  illustrates  that  moderate  waves  do  little 
mixing  and  serve  only  to  pump  the  thermocline  (and 
temperature  field)  up  and  down. 

5.  The  temperature  field  lags  the  wave  field  by  ap- 
proximately 180°;  but  the  lag  is  modified  by  as 
yet  undetermined  factors  and  depth.   The  spectral 
analysis  yielded  the  above  with  respect  to  surface 
waves  and  the  results  of  the  isotherm  plot  in 

Run  3  indicate  the  same  for  internal  waves. 

6.  The  model  demonstrates  the  high  probability  of 
statistical  independence  between  turbulent  and 
wave-induced  temperature  fluctuations. 

7.  The  turbulent  temperature  spectra  agrees  very  well 
with  inertial  subrange  theory  in  that  the  slopes  of 
the  log/log  spectra  are  very  near  -5/3.   That  the 
slope  is  nearer  to  -2  may  indicate  that  the  time 
series  was  non-stationary. 

8.  The  low  frequency  end  of  the  turbulent  spectra 
appears  to  have  different  slopes  depending  on 
whether  the  individual  run  had  internal  waves 
present . 
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quency and  wavenumber  spectra  were  calculated  and  a  correspondence 
between  the  spectra  was  noted  at  the  frequency  and  wavenumber  of  the 
surface  wave-induced  particle  displacements. 

The  Thornton,  Boston,  Whittemore  model  of  wave-induced  tempera- 
ture fluctuations  was  tested  and  found  to  model  the  temperature  spec- 
tra quite  well,  especially  in  a  narrow  band  of  frequencies  associated 
with  surface  waves.   The  turbulent  temperature  spectrum,  calculated 
as  the  difference  between  the  actual  and  wave-induced  spectra,  had  a 
slope  near  -5/3  above  0.1  Hz  and  more  negative  at  lower  frequencies. 
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